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Abstract — This paper introduces the Discrete Dithered Desyn- 
chronization (D 3 sync ) algorithm which is a decentralized Time 
Division Multiple Access (TDMA) technique in which a set of 
network nodes computes iteratively a conflict-free schedule so 
that each node obtains a portion of a frame that is an integer 
multiple of a fixed slot size. The algorithm is inspired by the 
dynamics of Pulse Coupled Oscillators (PCO), but unlike its 
predecessors that divide arbitrarily the frame among the nodes 
in the network, the D 3 sync allocates discrete resources among 
the network nodes. 

Our paper proves the convergence of the D 3 sync algorithm 
and gives an upperbound on the convergence time of the 
algorithm. 

Index Terms — Pulse coupled oscillators, desynchronization, 
and decentralized scheduling. 

I. Introduction 

Agreeing on a common timing (syncrhonization) and agree- 
ing on a Time Division Multiple Access (TDMA) schedule 
are related problems. In centralized scenarios there is a fixed 
infrastructure of nodes (masters) that optimizes the schedule 
and communicates them back to the other nodes (slaves). 
In infrastructure-less networks TDMA schedules are typically 
attained using decentralized graph coloring algorithms (e.g. Q 
and I3). In the first algorithm time synchronization between 
the nodes is presumed, and the second algorithm uses extra 
control messages to indicate the system state. 

Synchronization schemes themselves can be divided into 
these two categories. Centralized algorithms rely on the broad- 
cast signal coming from some fixed units, equipped with 
accurate timing devices; all nodes synchronize their clocks 
independently. The main example is provided by the Global 
Positioning System (GPS). In decentralized algorithms, in- 
stead, reference signals are generated at random or by a GPS 
signals, but then all users share information about the reference 
signals they sense, either by flooding the time stamp (e.g. Q) 
or by reaching consensus (e.g. l4l). 

Due to the issue of coverage, centralized algorithms have 
a considerable infrastructure cost. Decentralized methods are 
cheaper and less vulnerable to failures; however, they are 
slower and generally less accurate than centralized methods. 

In the context of infrastructure-less networks, great attention 
has been paid over the past ten years to biologically inspired 
synchronization algorithms based on the Pulse Coupled Oscil- 
lators models Q. These designs can overcome the necessity of 
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a reference signal for synchronization, since a reference signal 
emerges from the collective transmissions of the nodes. All 
PCO based algorithms, in fact, borrow a key implicit technique 
to communicate and update the scheduling decisions, which is 
through the emission of a special signal, called firing signal. 
The nodes in range of the firing node, update a local state 
variable as they detect a firing event. Another variant of time 
synchronization is desynchronization [6], which can be used to 
attain a TDMA schedule in continuous time in a decentralized 
fashion. More specifically in [6 ] the authors proposed the algo- 
rithm called Desync for a fully connected networks to attain 
a uniform TDMA schedule. Several different methodologies 
were utilized afterward to achieve deysnchronization, such as 
171 and (8). In the former, desynchronization is aimed via 
a graph coloring algorithm using a message-passing method 
with a constraint on the complexity of the messages. While in 
the latter, the proposed algorithm uses a PCO based method 
with a new kind of updates based on the firing signals received 
from all other nodes. Other variants of the desynchronization 
were proposed later with various conditions and objectives. 
In a Proportional Fair Scheduling (PFS) algorithm was 
proposed to achieve a time schedule in which nodes are 
assigned a fraction of the frame proportional to their relative 
demand. In fTOlL desynchronization is gained with another 
constraint on having a single node with a fixed time reference 
in the network. In general, compared to decentralized coloring, 
PCO based scheduling eliminates the need of synchronization 
to attain TDMA scheduling. However, most of the PCO 
based scheduling algorithms work under the condition of full 
connectivity in the network. In iTTTTl the authors tried to extend 
the Desync to multi-hop networks. 

The objective of this paper is to deal with the problem 
of providing a desynchronization algorithm that assigns not 
continuous but discrete portions of the frame. This modifica- 
tion is motivated by the fact that firing signals and modulated 
signals that are sent by the nodes are possibly of non-negligible 
duration compared to the frame size and, therefore, nodes 
cannot really use an arbitrarily short amount of the frame 
duration for their transmission. 

Our paper shows that, while quantizing the results of 
desynchronization leads to undesired fixed points for the al- 
gorithm, the algorithm we proposed, coined Discrete Dithered 
Desynchronization (D 3 SYNC ), converges almost surely to the 
desired resource allocation. A version of this protocol was 
proposed in our previous work in [12], but the effect of 



dithered quantization was assessed only by simulations. The 
analysis we carry out in this paper also provides bounds on 
the convergence time of the algorithm. Some of our analysis 
takes inspiration from the methods used in fT3l to prove the 
convergence of Quantized Consensus. 

This paper is organized as follows, in section [TTJ we intro- 
duce the D 3 SYNC algorithm. In order to analyze the charac- 
teristics of this algorithm, an equivalent consensus problem 



is defined in section III and it is proved that this problem 
converges to a fixed group of states. In section |IV| an upper 
bound for the worst case expected convergence time of the 
D 3 SYNC algorithm is computed. Simulation results are shown 
in section |VJ and a conclusion is made in section [VI| 

II. Discrete Dithered Desynchronization 
Algorithm 

Consider a fully connected network J\f with N nodes each 
with a discrete time counter for i = 1, . . . , N. Let the time 
unit for counters be A and assume that the node counters can 
only take values from the finite set 71 = {to, . . . , tl-i}, 
where Tj = jA. Like an alarm that is advancing towards 
its termination, the counters advance in time as ^(t) = 
A([t/A\ + x/ji (mod L)), where ^ £ {0, ...,L - 1} for 
all i = 1,2, ...,7V denotes their initial phases. Each node 
broadcasts a firing signal to the network upon overflow in its 
counter, i.e. when ^i(t) is reset to after LA periods. In the 
rest of the paper, similar to (6), the time that elapses between 
the firing of two consecutive counters is the time allocated to 
each node, in the order of firing. Without loss of generality 
we can relabel nodes to have < ipi < ip2 < • • • < ipN < L. 
It is implicitly assumed that no two counters can acquire the 
same value initially^ Nodes update their own counters when 
they sense firing signals from others and, in the absence of 
processing errors, desynchronization algorithms preserve the 
ordering in their evolution. Throughout the article we assume 
that all the additions and subtractions of the nodes indices are 
done modulo N, unless otherwise is stated. 

Because their initial position is random the nodes do not 
share a fair schedule, and the objective of D 3 SYNC update 
rule is to have the counters converge to a configuration 
for their counters within the frame approaching a uniform 
discrete TDM A schedule. It is conventional to describe the 
configuration of the counters over time around a circle that 
represents the frame duration, as shown in Fig. [T] Given 
that the order is based on the temporal proximity of the 
counters, we will refer to as time-neighbors nodes whose 
indexes differ by one unit. However, as can be noted from 
Fig. [T] the nodes geographical position is unrelated with the 
nodes relative counters position on the dial, which we assume 
to be completely random. For simplicity, we assume that all 
nodes can hear each other; strictly speaking, as clarified next, 
all we need is that nodes can hear their time-neighbors. 

! The assumption is that in the initial random configuration no two nodes 
have chosen the same counter. If this condition is violated there is a collision 
and the nodes will try again. 
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Fig. 1. Network Af with 6 nodes at t = t' , where (£') = 13A, ^ 2 (0 = 
14A, tf 3 (f) = 16A, ^ 4 (*') = 4A, ^ 5 (*') = 5 A, tf 6 (t') = 10A and 
L = 16. The grey nodes inside represents a random deployment of the nodes 
in A/", and the green nodes on the circle are their counters around a dial, 
which advance in time towards the finish line, at the north pole. 



The D 3 SYNC is a modification of the DESYNC algorithm. 
In Desync the update for node i occurs when node i — 1 
sends its firing signal at t = U-\, so that its counter is closer 
to the midpoint (in time) of the i — 1 and i + 1 counters: 

$i(tU) = + ^ (*i+i(*ti) + *i-i(*ti))> 

(1) 

where a G (0, 1), and denotes the time right after 

the updates resulted from node i — 1 firing signal. Because 
node i — 1 just fired, its counter is reset at time tf_ 1 , 
i.e. ^ i _i(t^L 1 ) = 0. Also, since only node i updates its 
counter in the network at U-\, the counter of node i + 1 
remains unchanged = \E^ + i(^_i). Therefore, ([T]) 

is equivalent to: 

= + ^-^*t+i(*i-i). (2) 

A simple way to jnodify this algorithm to obtain a discrete 
schedule so that ^(t^) G TL is to uniformly quantize the 
counter with the mapping Q(x) = min r . e 7^ \tj — x\. How- 
ever, as it will be discussed in Lemma [T] this approach may 
not converge to a desired fixed point. Hence, in D 3 SYNC we 
used a randomized Q, which is a dithered quantization lfT4l 
over Tl defined as follows: 



Q(x) 



mm Tj 

Tj€~L 



(x + v)\, 



(3) 



where x G K and v be a random variable, uniformly dis- 
tributed over (— A/2, A/2), v and x are statistically indepen- 
dent. Let Tj < x < Tj+i, then from equation Q distribution 
of Q(x) can be written as Pr{Q(x) = t^} 

- if i = 2, 

if z = .7 + 1, 
otherwise. 

(4) 
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The idea is analogous to the probabilistic quantization used 
in fT5ll to ensure the convergence of an average quantized 
consensus policy. 

Using dithered quantization ^ on (J^) ^(^L-J is mapped 
onto 72: 

*i(*ti) = Q(M4-i))- (5) 

Thus, Discrete Dithered De synchronization (D 3 SYNC ) algo- 
rithm works as shown in Fig. [2] and [3] 



Fig. 2. The D 3 SYNC algorithm for node i. When its counter overflows 
(Sending state), or when it receives the firing signal (Receiving state). Note 
that all the node needs is to follow the evolution of its own counter and wait 
until it hears a node fire after its own firing event. The use of indexes to 
refer to the nodes is for the convenience of our description, but it is totally 
irrelevant for the update. 

It should be clear at this point that, if there are no errors, 
the order of firings will be preserved. In fact, if initially all 
counters have at least a unit difference between each other, by 
using the semi linear update ^ the updated counter cannot 
cross the two fixed time-neighbors which are both on TL- 
In other words, the midpoint never falls inside an adjacent 
quantization bin of any two time-neighbors, and hence it is 
impossible to be quantized to any of their values. 

It is important to note that node i does not have direct 
access to other nodes' counters (i.e. and ^+i). The 

only information it acquires is the counter difference between 
its own counter and the node which last fired, prior to its own 
firing event, and hence immediately after it. In particular in 




Fig. 3. D 3 SYNC updates: Time evolution model of network Af with 6 nodes 
at t = t', where ^i(t') = 13A, ^ 2 {t') = 14A, tf 3 (*') = 16A, ^ 4 (V) = 
4 A, ^ 5 (*0 = 5 A, V 6 (t') = 10 A and L = 16. 



the D 3 SYNC algorithm, node i can compute the exact counter 
difference with node i— 1 and only an estimation of the counter 
difference with node i + 1. For the rest of this article, we 
assume that node i has exact information about the counter 
difference with the two time-adjacent neighbors (i.e. nodes 
i - 1 and i + if] 

The key ingredient of our incoming analysis, is a mapping 
of the counters used in D 3 SYNC onto a set of auxiliary 
variables. In fact, for the sake of the analysis it is preferable to 
consider the counter differences instead of the absolute counter 
values, since the sought TDM A schedule is equivalent to attain 
consensus on the counter differences. More specifically, let 
q(^) = (Qi(t)j • • • ? QN(t)) T denotes the difference between 
consecutive nodes counters, defined as follows: 

fc(*) = *»(f)-*i-i(t) (modL). (6) 

In order to rewrite the update equations based on the counter 
differences, after substituting ([6]) in ^ and ([5]): 

= Q(—2-qi(U-i) + -g-fc+ifa-i)). (7) 

It is clear from the definition in ^ and Fig. [T] that for any t 
the sum of all counter differences remain constant 

N 

J2 qi (t) = LA. (8) 

Hence, as qi and q i+ i are the only participants in the update 
equation ^ at time t = U-i, the sum of the two counter 
differences should be preserved after the update: 

Qi(4-i) + Qi+i^t-i) = Qi(U-i) + ft+i(*t-i). (9) 

As it was stated before, our final objective toward the 
implementation of a TDMA schedule over discrete resources 
is to divide a finite number of time slots equally among all 
the nodes in the network. But, this is only possible when 
L = £N for some integer £, and therefore each node will 
have £ slots (i.e. qi = £A for all i = 1, . . . , N). We define 
more precisely what we mean by a TDMA schedule in general 
in the following definition: 

Definition 1. TDM state: If L = £N + r for < r < N, then 
we denote TDM states as all subdivisions of the slots among 
the nodes in which r of nodes have i+1 slots and others have 
£ slots. 

In the following section we define an equivalent description 
of the evolution of the network state that eases the analysis 
of the behavior of the D 3 SYNC algorithm. The goal of our 
analysis is first to prove that the network dynamics will 
converge to one TDM state from any initial condition. We 
also provide a bound on the longest expected amount of time 
it takes to reach a TDM state, starting from the set of worst 
initial conditions. 

2 The assumption amounts to considering an ideal Physical (PHY) Layer, 
that detects and estimates perfecting the epoch of the firing signal. Simulations 
done in 1 9 1 have shown that the algorithm is robust and can recover from errors 
(false alarms or missed detection). 



D 3 SYNC Algorithm: 

Sending state: 

* Node i fires when its counter overflows (i.e. t = U 
such that = LA). 

* Node i resets its counter *&i(tf) = after firing. 
Receiving state: 

• Node i receives the first firing signal after its Sending 
state (from node i — 1) at £ = 

• *»(*ti) = G(a*»(*»-i) + 
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III. Analysis of the equivalent consensus problem 

As stated at the end of the previous section, in D 3 SYNC the 
values of the counters' differences in ^ at any firing event 
fully determine the evolution of the network Af schedule over 
time. It is assumed that all counters have a synchronized 
rise edge. Thus, qi (which is the counter difference between 
the pair of nodes i — 1 and i) may only change right after 
the firing events at U-\ or U-2, and it remains unchanged 
for all other firing events until the next round. Hence, we 
state that D 3 SYNC corresponds to the equivalent consensus 
problem described next. Suppose Q is a ring graph with N 
vertices V = ^2, . . . , v^}, each possessing a value qi for 
i = 1, 2, . . . , N. These vertices are connected to each other 
by the ring graph Q as depicted in Fig. |4| where a single 
edge, named as active edge, is distinguished from all other 
inactive edges. The two vertices connected to the active edge 
are called active nodes. The active edge also rotates counter 
clock- wise. Let each disposition of the active edge in time 
represents an event in Q denoted as {774,72,2,77,3,...}. For 
instance, if (^,^+1) is an active edge at rij, then 
will be an active edge at n J+ i. Now, if qi G Tl for all i and 
Y^f = iQi = LA, and also at each event active nodes perform an 
interaction on their values based on update equations similar to 
([7]) and ([9]), then this consensus problem would be equivalent 
to performing D 3 SYNC algorithm on network Af. In general 
if (vi,Vi+i) is an active edge at rij, then the stored values qi 
and q i+ i are updated as follows: 

, _/l + a 1 - a \ 

<6 = Q {-^^ + -^fc+V' (10) 



■Qi+i -Qn 



(11) 



where a G (0, 1), and q[ and q' i+1 represent the updated values 
of qi and q^i respectively and right before n J+ i. Note that for 
any x G Tl and y G M, it follows that Q(x + y) = x + Q(y). 



Thus, as qi+i G Tl, (10) can be written as 

1 



^ /l + a 



Qi+i 



2^)ft+i) 
fc+i))- 



(12) 



In the following parts, we investigate the behavior of the 
network's state q = (gi, . . . , gA/-) T in an evolving graph 
(5. Our objective is to show that the D 3 SYNC algorithms 
converges to the TDM states almost surely. The effect of 
uniform quantization in flQ| is discussed in the following 
lemma. 

Lemma 1. In the D 3 SYNC equivalent problem, if uniform 
quantization is used instead of the dithered quantization. Then 
there exists cases where the algorithm converges to a non- 
TDM state. 

Proof: The proof for this lemma is in Appendix [A] ■ 

Definition 2. Let be an active edge in graph Q at 

rij. Then, the possible interactions are: 




Active edge 



Fig. 4. The equivalent consensus problem for graph Q with N nodes. The 
active edge is currently (vi, fj+i), and it will rotate in a counter clock wise 
direction at each time step. 



1) Null: When q[ = qi and q[ +1 = q i+1 . 

2) Swap: When q[ = q i+ i and q' i+1 = qi. 

3) Compression: When \q[ — q' i+1 \ < \q% — 

The following lemma discusses the behavior of the interac- 
tions in the D 3 SYNC equivalent consensus problem. 

Lemma 2. Let (1^,1^+1) be an active edge in graph Q and 
^ T L , then: 

1) If qi = qi+i, with probability one q^ = q' i+1 . 

2) If \qi — qi+i\ = A, the updated values will have the 
same difference afterward \q[ — q' i+1 \ = A. 

3) If \qi-q i+ i\ > A: 

a) If qi > q i+1 : 



q'i<qi and q' i+1 >q i+ i. 



(13) 



b) If qi < q i+1 : 



q'i>qi and q' i+1 < q i+1 . (14) 

4) After each interaction \q[ — q' i+1 \ < \q% — q%+i\. 

Proof: The proof for this lemma is in Appendix [B] ■ 
This lemma characterizes all the possible interactions in 
Q. Let (vi,Vi+i) be an active edge in Q. If both of the 
active vertices have the same values, with probability one the 
null interaction will happen. Also, if the values of the active 
vertices have more than one unit difference (like in Lemma [2] 
case 3), the smaller one will increase and the bigger one will 
decrease, i.e. the interaction is a compression. It is important to 
note that in general the difference between the values would 
decrease with a positive probability, which depends on the 
choice of a. 

Corollary 1. If active nodes have a difference of A in their 
values (Lemma [2] case 2) two cases may happen: 1) a null 
interaction, with probability or 2) a swap interaction, 

with probability 

Based on the results of Lemma [2| it is simple to see that the 
properties of the D 3 SYNC equivalent consensus problem are 
similar to the Quantized Consensus in (T2l on the ring graph Q. 
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Reaching an average quantized consensus over Q is equivalent 
to having the desired discrete desynchronization in the original 
network M . There are two important differences, however, 
between the D 3 SYNC algorithm and the Quantized Consensus: 

1) The edge selection in our algorithm is deterministic (i.e. the 
active edge deterministically rotate counter clockwise), while 
in the Quantized Consensus the edges are chosen at random. 

2) In our algorithm the updates are randomized, similar to 
[Q2L while in Quantized Consensus there is a deterministic 
operation over the chosen edge. The random evolution of the 
network state q = (gi, . . . , Qn) T in Q is a Markov chain. 
In fact, to specify the distribution of the next state q' all is 
needed is the current position of the active edge and the current 
network state q (in particular the values of the active nodes). 
Also in the equivalent consensus problem YliLi Qi ~ ^A 
and qi > A (based on the constraints of the D 3 SYNC al- 
gorithm). Consequently, there are finite number of possible 
states for network Q. Because of the deterministic movement 
of the active edge, D 3 SYNC equivalent consensus problem is 
a cyclo- stationary Markov process. In contrast, the Quantized 
Consensus in fT3ll and the average consensus method in fT5l 
form stationary Markov processes, and hence the evolution can 
be modeled as a homogenous Markov chain. This difference 
is particularly important in characterizing the convergence 
time, which is more complex for D 3 SYNC than for Quantized 
Consensus. Establishing almost sure convergence is simpler, 
and done to prove the following results. 

Theorem 1. In the D 3 SYNC equivalent consensus problem on 
network Q, TDM states are absorbing states. 

Proof: Let S L = {q : |^ - qj\ < A, Mv u Vj G G, ||q||i = 
LA} denotes the set of all TDM states in Q. It suffices to 
show that once the network is in a TDM state (i.e. q G Sl), 
the updated network will also be in one of the TDM state (i.e. 
q' G Sl)- Suppose q G Sl and (1^,1^+1) is the active edge. 
Then for those interactions in which both vertices of the active 
edge have equal value (i.e. qi = (ft+i), based on Lemma [2] case 
1> Qi = Qi+i- Therefore, q' = q G Sl- However, for those 
interactions in which \q^ — = A, based on Lemma [2] 

case 2, \q[ — q f i+1 \ = A, which means either a null or swap 
interaction will occur for and qi-\. In either case, all the 
other inactive vertices will have their values unchanged and, 
hence, q' G Sl, which proves the statement. ■ 

Theorem 2. In D 3 SYNC equivalent consensus problem on 
network Q, the TDM states are the only absorbing states and 
the network converges to the TDM states almost surely. 

Proof: Let q(no) be the initial state of the network Q. 
We want to show that after sufficient number of interactions, 
with probability one, the network reaches the absorbing state: 

lim Pr{q(n) G S L } = 1- (15) 

n— )*oo 

Let Rg(n) = max^(n) — mirij qj(n) denote the range 
between the maximum and minimum qi (n) in network Q at the 
nth time step. Based on Lemma [2] case 4, after each interaction 



the difference between the values of the two active vertices 
is a non-increasing quantity, i.e. Rg(n + 1) < Rg(n). In 
particular, Rg(n + 1) < Rg(n) when there is a single node 
with the maximum (or minimum) value in the network which 
participate in a compression. Also, < iiV( n ) ^ A only 
when the network is in Sl', in fact, Rg(n) = only when 
L = £N and Rg(n) = A otherwise. Instead, Rj\r(n) > A 
for all the states which are not in Sl- Therefore, to prove 
the theorem, it is sufficient to show that while q(n) ^ Sl, 
Pv{Rg(n + k) < Rg(n)} > 0, for some k. In other words, 
it is possible for the network to move into a state with lower 
range after finite number of interactions. Suppose the network 
state is not in Sl- As the active edge rotates over the entire 
network from this initial state, there exists a time step n* 
in which one of the two vertices of the active edge has the 
value q max = max^(n*). We denote that active vertex as 
vm- Since, q(n*) ^ Sl, there should exist at least one other 
node in Q with value lower than q max — A, otherwise the 
initial network state considered must be in Sl, which is a 
contradiction. Now moving counter clockwise from vm, let 
Vj be the closest node such that: 

j = M - argmin{g M -/c < ^max - A} (mod N). (16) 

k 

Then all the nodes between vm and Vj (moving counter 
clockwise) have values equal to either q max or q max — A. Based 
on Lemmas [2]- 1 and [2]- 2, with positive probability the value 
of vm can swap until qj+i = g max , and then, with positive 
probability, a compression will happen on the active edge 
(vj,Vj+i). Depending on the number of nodes having the 
value equal to g max , by following the previous steps, one can 
prove that, with positive probability, the number of nodes 
with maximum value is reduced by one, until all of them 
have participated in a compression, which is the case where 
Rj\T is strictly decreased. Thus, independent of the initial 
state, with positive probability Rjy is decreased step by step 
and with a finite number of interactions, until it reaches the 
minimum value which correspond to an absorbing state. Since, 
the sum of the qi is preserved during each interaction, it is not 
necessary to keep track of the minimum values, but the same 
process happens on the minimum values too. 

It was shown that, starting from any state that Rg > A (i.e. 
states which are not in Sl), with positive probability Rg is 
strictly decreased after sufficient number of interactions. Thus, 
the only possible absorbing states are Sl which < Rg < A. 

■ 

IV. Analysis of the Convergence Time for the 
D 3 sync Algorithm 

Up to this point, it was shown that irrespective of the 
initial network state (i.e. q(no)), after a sufficient number 
of interaction the system reaches a TDM state in Sl almost 
surely. It is now of interest to evaluate the expected number 
of interactions needed for a system to reach a state in Sl- 

Let S = {q : Zf=iQi = LA, Qi > A,W, G V} 
be the set of all possible states in the network Q, and 
V(q) = ||q-^l|| 2 /A 2 , where q = ± Qi . A similar 
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function is defined in fT3l to characterize the dynamics of the 
average consensus with quantized values. Similar to Lemmas 
4 and 5 in fT3lL it can be shown that: 

Lemma 3. If q G <S, V(q) has the following properties for 
the network Q with ||q||i = LA: 

1) max qG5 V(q) < (L — N) 2 N/4. 

2) After each interaction: 



V(q') = V(q) if \ qi -q i+1 \<A, 
V(q / ) < V(q) - 2 otherwise. 



(17) 



Proof: As E^Ii^W = LA > and Qi(t) > A > 
max^ qi(t) < (L — N + 1)A, which can only occur when all 
other values are equal to A. Thus, the first property is achieved 
by modifying Lemma 4 in fT3l by the following bounds. The 
second property is another representation of Lemma 5 in fT3l 
for this scenario. ■ 
V(q_(n)) is a positive function which is decreasing in n, 
and will reach its minimum once the system reaches the TDM 
states (i.e. q(n) G Sl)- The intuition behind this fact is that 
once the system is in Sl, all qi(n) are either equal or can only 
have a unit slot difference, then from the second property of 
Lemma[3]F(q) does not change anymore. Also, from Theorem 
[2] the system will eventually reach the TDM states in Sl, 
equivalently V(q) will reach its minimum. 

As it was mentioned earlier, the dynamics of the network 
Q can be modeled as a Markov chain. Let T(no) = min{n > 
no : q(no) ^ <S/,,q(n) G Sl} be the absorption time 
of the Markov chain. Let Ti(q, n*) be a random variable 
defined as the number of interactions, starting from step n* 
and state q(n*), until the first compression in the network 
Q. Let f(no) = max qG s E{Ti(q, no)} be the maximum 
expected number of interactions one has to wait until network 
Q experiences a compression. Then based on the results from 
Lemma |3j it can be shown just like in lfl~3ll that: 

max E{T(n )} < f(n )(L - N) 2 N/8. (18) 

q(n )e*S 

The intuition behind the inequalit yjT8] ) is as follows: Based 
on the second property of Lemma [3]depending on the initial 
state q(no) ^ Sl, for n > no, V(q(n)) will decrease only if a 
successful compression happens. In the worst case, V(q_(no)) 
is maximum and at each compression it is decreased by the 
smallest possible amount, which is by 2. Also in this scenario, 
two consecutive compressions happen after at most T(n) 
number of interactions (in expectation), given the fact that 
at the nth time step a compression occurred. Next we provide 
the main analytical result of this paper: 

Theorem 3. In the network Q with TV nodes and update 
parameter a G (0, 1), the worst case maximum expected 
number of null and swap interactions until a compression 
occurs, T (no), can be calculated for large N as: 



f(n ) = o(g(a)N 3 ), 



(19) 
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Fig. 5. Convergence of the D 3 SYNC algorithm over two scenarios for a 
network J\f with 6 nodes and a = 0.2: (a) L = 60, where there exists a 
unique absorbing state, (b) L = 57, where \Sl\ = 20. Both scenarios are 
converged after 6 rounds of firings. 



Proof: The proof of this theorem is in Appendix |C] ■ 
Based on Theorems [3] and (18), the greatest expected 
number of interactions to reach a TDM state is upper bounded 
by 0(g(a)(L - N) 2 N 4 /S). Also, g(a) is a constant factor 
which is an strictly increasing function of a for a G (0, 1). As 
it will be shown in simulation results in Section |V| in general 
this upper bound is not tight for the D 3 SYNC algorithm, and 
it represents an upper bound for the worst case possible. 



where g(a) = 2 4{i-a) • ^ e exact vame of T{uq) for any N 
can be found in equation (41) in Appendix [C] 



V. Simulation Results 

In this section performance of the D 3 SYNC algorithm in 
TDMA scheduling and also its convergence time is assessed 
by using computer simulations. 

The D 3 SYNC algorithm has a single parameter a to be 
chosen. It represents the algorithm's inertia in updating the 
current phase counters based on the received firing signal. As 
it is clear from ([19]), large a shows more resistance toward 
the update, and hence its convergence time is greater in 
comparison with a smaller a. 

In Fig. [5] convergence of the D 3 SYNC algorithm to the TDM 
states is shown under two different scenarios. Network Af is 
consisted of 6 nodes with a = 0.2. In (a) L = 60 and the 
network has a single TDM state (i.e. = 1)> while on (b) 
L = 57 and hence there is not a single TDM state. After six 
round of firings, where each round is defined as a complete 
cycle of firing events in the network and by all of the nodes, the 
network is absorbed into Sl- There are \Sl\ = 20 different 
TDM states and the network changes its state among these 
states afterward. 

Number of interactions until absorption for the D 3 SYNC al- 
gorithm is depicted in Fig. [6]-a. As it is clear from the figure, 
the upper bound in ([18]) is very pessimistic, because it is 
based on the performance of the worst possible case as for the 
initial state. Theorem [3] found an expression for the expected 
number of interaction until the next compression in network 
Q. As it is explained in Appendix [C] this is the case where 
all nodes have a common value, except for the two of them, 
one of which a unit higher and the other a unit lower than 
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Fig. 6. (a) Network J\f with N = 10. Number of interactions until absorption 
versus L averaged over 5000 different random initial states, (b) Network 
Af, with all nodes having a common value, except for two of them, one 
of which a unit higher and the other a unit lower than the common value. 
Number of interactions until absorption versus N is depicted averaged over all 
possible initial positions and each with 250 independent trials. The dashed 
lines represent the theoretical number of needed interactions computed in 
Theorem [3] 



the common value (also shown in ([30])). The average number 
of interactions until absorption for this case is shown in Fig. 
[6]-b. The evaluated absorption time in Appendix [C] which was 
discussed in Theorem [3] is also shown in this figure by dashed 
lines. 

VI. Conclusion 

In this paper we showed that by using the D 3 SYNC algo- 
rithm, desynchronization and , hence, TDMA scheduling is 
possible over discrete resources and in a decentralized manner. 
This approach can also be used to reach TDMA in more 
realistic situations, such as allocating resources based on each 
node's demand, or also considering the case where rise edge 
times of counters are not matched. 

Appendix A 
Proof of Lemma[TJ 

It was stated in Section [TT] that by using uniform quantization 
on ([2]), the algorithm may not converge. We continue our 
discussion on the D 3 SYNC equivalent consensus problem. Let 
qi = q i+ i + A, then by replacing the uniform quantization in 
( fT2] ) (note that Q(-) denotes the uniform quantization): 

q'i = Qi+i + Q(^^A) = q i+1 + A. (20) 



where the right hand side equality is because for any 
(0, 1), \ < ±±« < 1. Similarly, if q { = q i+1 - A, then 

1 



q[ = qi+i + Q( A) = q i+1 - A. 

In a network with N = 4, q = (A, A + 1, A + 2, A • 
a fixed point, but it is not a TDM state. 



(21) 
If 



is 



Appendix B 
Proof of Lemma[2] 

To prove statement 1) is true we observe that, since 
Qi,Qi+i £ Tl, by replacing q i+1 = q { in equation dTOb, 



q'i = Q(Qi) = Qi, an( i consequently from (11) q f i+1 = q i+1 . 

To prove statement 2), we consider equation ( [T2| ). Two cases 
may happen. First, suppose qi — q^i = A. Then q[ = q i+ i + 
q(^a), and since < a < 1, then \ < ^ < 1. Thus: 



2 

1+q 
2 * 



qi+i with probability 

q i+ i + A with probability 

In either case, it can be inferred from ( [TT] ) that: 

, = {qi tfq , i = q i + 1 , 

q%+1 \qi-A if g( = + A. 

Thus in general, 



■Qi+i 



with probability ^f- , 

l+OL 

2 * 



with probability 



(22) 



(23) 



(24) 



Second, assume q i+1 — ^ = A. Note that Q(—x) 
because of the fact that the distribution only depends on the 
distances and hence any transformation on x which keeps 
distances unchanged (like taking the complement) has the 



same exact effect on Q(x). Using this fact, and equation ( 12), 
the update can be written as q[ = qi+\ — Q^^^A^. With 
similar reasoning as before, the updated difference can be 
finally written as: 



■Qi+i 



with probability 
with probability 



l-q 
2 ' 

1+a 
2 * 



(25) 



Therefore this proves that, if \qi — q^i | = A, after the update 
W i -q' i+1 \ = A. 

For the proof of statement 3), suppose without loss of 
generality that \qi — qu+i | = mA, m > 2 and qi > qi+i. 
By using equation (10), q[ = Q^-^qi 
and hence, 



mA)j, 



qi 



a 



mA 



(26) 



Now, since ^mA is a positive value, the right hand side is 
greater than or equal to the left hand side, an hence q[ < q%. 
The sum of the two values updated will be preserved after the 
update, thus, from equation ( [TT] ) it can be readily concluded 
that q' i+1 > qi+i and the first part of the statement is proved. 
For the second part, the proof is the same as the first one 
because of the symmetry between the two cases. 

Statement 4 combines the previous ones. If qi = q i+1 or 
\qi — qi+i\ — A. If \qi — _^_mA, m > 2, suppose 

i\ and 



G that qi > qi+\. Thus in equation (26), if m! = [- 



V 



(1— a)m 



m' the second term in the right hand side can 



be written as 



-rriAj = 



ra'A 



with probability 1 — p 7 , 



(m' + 1)A with probability p f . 



(27) 
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Thus, equation ([26]) can be manipulated as follows: 



N-l N-2 



i + 1 



Qi ~ Qi 



-m'A with probability 1 — p' 

~(m' + 1)A with probability p' . 



(28) 



A linear combination of equation (11) and ( 28 ) is 



Qi ~ Qi+i 



2m') A with probability 1 — p f , 
2m' — 2) A with probability p' . 

(29) 
< mA 



Since < ^ < \ and < ra' < f , |^ - q' i+1 \ 
qi — qi+i. Similarly, if qi < q i+ \ by the same approach it 
can be shown that \q[ — < mA = q i+1 — qi, and this 

completes the proof. 

Appendix C 
Proof of theorem[3] 

For the proof of [3] we imitate the clever methodology used 
in fT3l to calculate T(n$). By definition T (no) is the expected 
number of interactions needed, in the worst case scenario, for 
network Q to have its next compression, starting from iteration 
no. For a compression to happen, it is necessary that the 
active edge (1^,1^+1) is such that \qi — <ft+i| > A. Thus, 
the worst case (in expectation) is when there are the least 
possible number of potential nodes in the network Q that can 
participate in a compression, i.e. only two nodes V{ and Vj 
have values \qi — qj\ > A and all remaining values are equal. 
In this case, to fall in the absorbing state, one would have to 
wait for a compression of the two outliers, while in all other 
interactions the value of V(q) will be unchanged. It is clear 
that any other initial state increases the probability of having 
a compression, and therefore this is the worst case scenario. It 
should be noted that if J2^ =1 qk(^o) — LA, the stated worst 
case scenario happens only when L = £N and all other cases 
(i.e. L — £N + r, 1 < r < N) require less time in expectation. 

Let us now focus on the worst case scenario, where there 
exists two nodes vi and Vj such that: 



(^ + 1)A k = i, 
(£-l)A k = j, 
£A k^i,j. 



(30) 



Any compression of the values that differ by 2 A from each 
other (initially at Vi and Vj) will lead to a TDM state, which is 
an absorbing state. Otherwise, the outliers can change position, 
but there will always be only two of them. Thus, the evolution 
of the network can be viewed as a Markov chain with a 
single absorbing state (q = £A1) and N(N — 1) transient 
states that are all possible permutations of positions for the 
outliers around the ring graph Q. In the following we study 
this Markov chain, and derive how long it takes on average for 
these two values, which are outliers, to get into a compression 
thereby settling the network in a TDM state. The state space 
is defined as S = £ U T, where £ = {(m, n) : 1 < m < 
N — 1,0 < n < TV — 1}, and m represents the relative 
position of the node with value {£ + 1)A with respect to the 
node with value {£ — 1) A (e.g. i — j (mod N) in the case of 



1 ^— m 

© r 




Fig. 7. The Markov chain for analysis of the absorption time of the proposed 
problem. The state space is S = SUJ 7 . Transition probabilities are also shown 
on each arc. The expected number of interactions starting from each state is 
written inside the corresponding circle. 



(pO])), and n represents all the N possible cases of the current 
active edge, where if qi = {£ — 1)A then the active edge 
(vi-i,Vi) is assumed to be labeled as and others are labeled 
until N — 1. T denotes the absorbing state. So, without loss 
of generality we have fixed our reference to the node with 
value (£ — 1)A and the current state is identified based on 
that. In Fig. [7] this Markov chain is shown, and the circle in 
row n and column m, represents the state (m, n). In Fig. |7j 
the expected number of interactions until absorption starting 
from each state, is written in the corresponding circle. The 
transition probabilities of the Markov chain can be derived 
based on the update equations ( fT0| ) and fTT] ). Recall that there 
is only one possible compression which can occur on the two 
outlier values, and it is possible only if both of their vertices 
are connected to the active edge. Based on the outlier values 



from (30) and replacing them into equations (10) and (11 
the probability of a compression in this case is 1 — a and, 
hence, the probability of a null interaction is a. If the active 
edge has two nodes in consensus, the state moves forward 
with probability one to activating a new edge with a distance 
between the outliers unchanged. For all other cases where the 
interaction is between an outlier value and non-outlier one, we 
know from Corollary [l] the probability of a swap is p = 
and the probability of a null interaction is 1 — p. There exists a 
symmetry in this Markov chain, which can be easily assessed 
by the choice of the shown variables x,y and z and the path 
through the absorbing states T in Fig. [7] The expected number 
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TV 4 (a + l) 2 + N 3 (a + l) 2 + 127V 2 (a - l) 2 - 247V(a - l)(2a - 1) + 24(a - 1) 

247V(1 -a)(l + a) 



(41) 



of interactions until absorption can be calculated as: 



Also, by doing some algebraic manipulations we have: 



N-3 
i=0 



N 



2 

J=0 



7V-4 N-3-k 



£ £ 4 WW- 1 )) (3D 



/c=0 £=1 



Based on the Markovian property, z| = 1 + and z£ = 

1 + for all valid fc and £ By eliminating all the z variables, 
the following equations can be written for all x and ?/: 



2/o - ayN-2 1, 

Xfe - (1 - p)2/ fc - px N - k _ 3 = 1 + fcp, 

2/fc' - PZ/fc'-l - (1 - P)XN-k'-2 =k' - 



(A/ - l)p, 



(32) 
(33) 
(34) 



for all < fc < TV - 3 and 1 < fc' < N - 2. Summing over 
all k and fc' in ( [33] ) and ([34]), together with ( [32] ), and 2/jv_2 
can be calculated as: 



a(7V-2)(7V + l) + 2 N(N-1) 

yo = ^ v ,vn-2 - 



2(1 -a) 



2(1 



a) 



(35) 



The other unknowns can be calculated recursively for 1 < k < 

N/2 as: 



+ (Vk-i +2/j\r-fc-i)/2, 



_ 2(l-p) 

x N -k-2 = (ajfe-i - Ffe - (1 - p)yk-i)/p, 

Vk = (k - (k - l)p) + py/fe-i + (1 - p)x N - k -2, 

VN-k-2 = (yN-k+i ~ (1 - p)^/c-i - Gk)/p, 



(36) 



where D k = (N - k - l)(2p - 1) + F fc = 1 + (A; - l)p, 
and Gk = (N — k — 1) — (N — k — 2) p. xo can be calculated 
based on ([36]) and by using y and 2/at-2 as: 



x 



A^ 2 (a + l) 2 + A/"(3a 2 - 6a - 1) - 10a 2 + 4a + 6 



4(l + a)(l-a) 



(37) 



Considering all of the solutions by equations d36j, it is not 
hard to see that x k - x k . x = ^Hyito-^ , 
and y k - = ^1)^2^)^-1))-^ ^ by 

using xo and yo we can calculate ^2 i=0 x % and Ylj=o % as 
follows: 

N-3 

J2 Xi = (JV - 2)x (38) 

i=0 

(TV - 2)(7V - 3)(N(N + l)(q + 1) + 6(1 - a)) 
+ 24(1 - a) ' 



AT-4 N-3-k 



N-4N-3-k 



/e=0 
AT-3 



/c=0 ^=1 



fc(fe+l) 



+ (7V-2-fe)x fc _ 1 



(40) 



Substituting ( [38] ), ( [39] ) and ( |40| ) into pT] ), the expected number 
of interactions until absorption (T) can be computed as in 
equation (41) at the top of the page. It can also be seen that 
for large N, f = 0(g(a)N 3 ) where g(a) 



24(1- a) ' 



N-2 

J2yj = (N-l)y (39) U5] 

3=0 

(N - 1)(N - 2)(7V 2 (a + 1) - 3A^(3a - 1) + 6(1 - a)) 
+ 24(1 - a) " 
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